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Abstract 



We investigate possible signatures of black hole events at the LHC in the hypothesis that 
such objects will not evaporate completely, but leave a stable remnant. For the purpose of 
defining a reference scenario, we have employed the publicly available Monte Carlo generator 
CHARYBDIS2, in which the remnant's behavior is mostly determined by kinematic constraints 
QQ and conservation of some quantum numbers, such as the baryon number. Our findings show 

that electrically neutral remnants are highly favored and a significantly larger amount of missing 
transverse momentum is to be expected with respect to the case of complete decay. 

rn 

O 1 Introduction 

• • Models with large extra dimensions [H [2] allow for the fundamental scale of gravity Mq to be 

as low as the electro- weak scale {Mq ~ ITeV), and microscopic black holes (BH) may therefore 
^ be produced in our accelerators (for reviews, see, e.g., Refs. [3j). Once a BH is formed, and the 

balding phase is over, the Hawking radiation |4| sets off. The standard description of this effect is 
based on the canonical Planckian distribution for the emitted particles, which implies the life-time 
of microscopic BHs is very short, of the order of 10~^^ s [5]. This picture (restricted to the ADD 
scenario ) has been implemented in several numerical codes [6l [TJ El EJ [IHl HH |12] , which have 
been mainly designed to extract information that will allow us to identify BH events at the Large 
Hadron Collider (LHC) (see Ref. [T3] for a search of BH signatures with the first data collected by 
the LHC and Ref. jl4j for some theoretical criticisms on the validity of that analysis). 

It is important to recall that the end-stage of the BH evaporation remains an open issue (see, 
e.g., Refs. [IS Ull ITT] ) , because we do not yet have a confirmed theory of quantum gravity. In fact, 
the semiclassical Hawking temperature grows without bound, as the BH mass decreases, which can 
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be viewed as a sign of the lack of predictability of perturbative approaches This is an important 
issue also on a purely experimental side, since deviations from the Hawking law for small BH mass 
(near the fundamental scale Mq) could actually lead to detectable signatures. 

Our ignorance of the late stages of the BH evaporation is usually bypassed in the numerical 
codes by letting the microscopic BH decay into a few standard model particles as an (arbitrary) 
lower mass is reached, but the possibility of ending the evaporation by leaving a stable remnant 
has also been considered [201 ED El] • The purpose of this report is precisely to study the case in 
which BHs stop emitting and leave a long-lived (or even stable) remnant of mass M ~ Mq. 

For this purpose, we have employed the Monte Carlo code CHARYBDIS2 [T^], which is presently 
the only generator which natively supports remnant BHs, as we shall describe in Section [3j In 
Section |4j we report results from test runs and describe in detail both the primary emission and 
final output obtained using Pythia [23]. We shall use units with c = h = 1 and the Boltzmann 
constant /cb = 1- 



2 Black hole production and evolution 

In the ADD scenario [Ij , a microscopic BH with horizon radius shorter than the typical size of the 
extra dimensions is usually approximated by a (4 + (i)-dimensional solution of the vacuum Einstein 
equations. In the simplest case of a non-rotating BH, this leads to the following expression for the 
horizon radius, 

1 /M\^/ 8r(^) \^ 



where M is the BH mass, T the Gamma function, and the temperature associated with the horizon 
is given by 

At the LHC, a BH could form by colliding two protons. The total BH cross section can be 
estimated from the geometrical hoop conjecture [23] as 

cj(M) 7ri?H . (2.3) 

In order to determine the total production cross section, this expression must be rescaled according 
to the parton luminosity approach as 



da 
dM 



= ^H(^b^BH;s = M^) , (2.4) 

pp-^BH+X 



where a and b represent the partons which form the BH, is their centre-mass energy and 

E 

a,b 



dM S ^ JaP/s \SXa ' 



^For reviews of the more consistent microcanonical description of BH evaporation in which energy conservation is 
granted, see Refs. [181 119| . 
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where fi{xi) are parton distribution functions (PDF) and ^/s the LHC centre-mass collision energy 
(up to 7TeV presently, with a planned maximum of 14TeV). Of course, semiclassical BHs will only 
form above a minimum mass (presumably) larger than Mq, [25] . 

Once formed, the BH begins to evolve. In the standard picture, the decay process can be 
divided into three characteristic stages: 

Balding phase: the BH radiates away the multipole moments inherited from the initial config- 
uration and reaches a hairless state. A fraction of the initial mass will also be radiated as 
gravitational radiation. 

Evaporation phase: the BH loses mass via the Hawking effect. It first spins down by emitting the 
initial angular momentum, after which it proceeds with the emission of thermally distributed 
quanta. The radiation spectrum contains all the Standard Model (SM) particles, (emitted on 
our brane), as well as gravitons (also emitted into the extra dimensions). For this tags, it is 
crucial to have a good estimate of the grey-body factors [261 127j . 

Planck phase: the BH has reached a mass close to the effective Planck scale Mq and falls into the 
regime of quantum gravity. It is generally assumed that the BH will either completely decay 
into SM particles [S] or a (meta-) stable remnant is left, which carries away the remaining 
energy [20j . 

We will here focus on the possibility that the third phase ends by leaving a stable remnant. 

3 Remnants and a Monte Carlo code 

Several Monte Carlo codes which simulate the production and decay of micro-BHs are now available 
(see, e.g. Refs. [H [71 [51 [U [lOl [TTl [12] ) . The only code which presently implements the description of 
a remnant is CHARYBDIS2 [T2] Such a description is mostly based on kinematic and quantum 
number constraints, so that the remnant will necessarily have zero baryon number. The production 
and evaporation phases are better modeled, from a theoretical point of view, although they could 
still be improved. For example, the production phase is described using bounds from trapped 
surface formation, which might lead to overestimate the amount of missing energy in gravitational 
radiation [28j. Further refinements have also been investigated for the Hawking spectra [29^ llOj. 

A list of the relevant adjustable parameters in the code and the values we used in the simulations, 
is given in Table [T] In the study we present here, we simulated the different scenarios in pp collisions 
with total centre-mass energy of both y/s = 7 TeV |^ and 14 TeV. 

It is particularly important to describe in details the role of one of the kinematic constraints, 
namely KINCUT. When set to TRUE, the generator will not allow the emission of Hawking quanta 
that reduce the BH mass M below Mq, whereas KINCUT=FALSE will allow the last decay particle 
to have an energy as large as the BH's (but no more). Let us then see what will likely happen 
for the last emission, when the BH has already evaporated down to a mass Mle ^ Mq. The 
generator will try to simulate an Hawking particle with energy cj determined by the Planckian 
distribution with a temperature Th = Th(Mle) and small enough to ensure energy conservation. 
Roughly speaking (that is, neglecting the BH kinetic energy), this means Th(Mle) ^ uj < Mle- 

^This code has been used to study the remnants of non-commutative BHs in Ref. [55] 

^This case was mostly included for a comparison with the CMS analysis in Ref. [13], although there are arguments 
which show that 7 TeV might not be theoretically enough to produce semiclassical BHs [14| . 
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Parameter 


Description 


Used values 


Mg 


Fundamental gravitational scale 


2-4 TeV 


D 


Total number of dimensions 


6, 8 


MINMSS 


minimum initial BH mass 


2 X Mg 


MAXMSS 


maximum initial BH mass 


Vs 


RMSTAB 


main switch for stable BH remnant 


TRUE/FALSE 


KINCUT 


kinematic condition for BH remnant 


TRUE/FALSE 



Table 1: Most relevant parameters of CHARYBDIS2 and their values for the different simulated 
samples used in this study. 



Since the temperature (2.2) increases for decreasing M, it is likely the attempted last oj will carry 
away enough energy to reduce M below Mq. If KINCUT=TRUE, this process is forbidden and 
the generator simply stops and leaves a remnant of mass Mle > -^g- If KINCUT=FALSE, the 
process is allowed and the remnant mass will be M ~ Mle — < Mg- The effect of KINCUT 
is clearly visible in Fig. [T| which shows the remnant mass distributions for samples generated at 
^ = 7 TeV with Mq = 2 TeV, in L» = 6 and 8. If KINCUT=FALSE the last decay step, which 
allows the remnant to access masses lower than Mq, typically releases a lot of energy, producing 
low mass remnants. On the other hand KINCUT=TRUE forbids remnant masses below Mg, thus 
producing an increase of events as the remnant mass approaches Mg and long tails towards higher 
values. These scenarios are in agreement with the results we shall show in the next Section. 

From a theoretical point of view, one expects a remnant mass around Mq (see, e.g., Ref. }25|). 
The above scenarios are therefore, in a sense, complementary, with the KINCUT=FALSE case 
perhaps less sound. 



4 Simulations 

The standard setting for CHARYBDIS2 simulates the evaporation of a microscopic BH according 
to Hawking's law, until the black hole mass M reaches the minimum value, which is assumed 
equal to the fundamental scale of gravity Mg- The BH subsequently decays into a few bodies 
(A^f = 2,..., 6) simply according to phase space. As described in the previous Section, we are 
instead interested in the study of scenarios with a residual stable BH remnant in the final state. 

In order to compare the output of the different scenarios (no stable remnant, stable remnant 
with KINCUT=TRUE/FALSE) we have generated Monte Carlo samples varying the parameters 
listed in Table [T] and assuming two different centre mass energy for the LHC, namely 7 and 14 TeV. 
As pointed out in Table [l} the mass of the produced BH has been required to be at least 2 Mq in 
order to avoid the region where the classical theory, used to evaluate the production cross section, is 
no longer valid (but see also Ref. [14J). The other parameters have been fixed to the CHARYBDIS2 
default values. Table [2] summarizes the characteristics of the generated samples, including the tt 
production process which is one of the most relevant backgrounds. The table also contains the 
event yields normalized to an integrated luminosity of Ifb^^ for few basic selection requirements 
applied at the hadronization level, as described in details in Section [5} 
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Figure 1: Remnant mass for different values of the KINCUT parameters. 



4.1 Primary Emission 

Let us first examine the "primary" BH emission, namely the particles produced by direct BH 
evaporation, before parton evolution and hadronization are taken into account. This emission is 
the direct output of CHARYBDIS2. Fig.[2]shows the produced final state particles and the particle 
multiplicity for y/s = 7 TeV, Mq = 2 TeV, in D = 8 and for the different remnant conditions. The 
final state particles, listed in Fig. [2| panels a), c) and e), are grouped in SM particles (the first 
three blocks contain quarks, leptons and bosons, respectively), gravitons (G), neutral (R'^) and 
charged (R^) BH remnants. The remnants are neutral in the large majority of the cases, with only 
a few percent of the events ending in a charged remnant. The neutral remnants behave as weakly 
interacting massive particles (WIMPs), producing a striking signal of missing transverse energy 
and a topology quite different from the usual case of total BH evaporation. 



4.2 Final Output 

The primary output includes unstable particles, along with quarks and gluons that have to be 
processed by a parton shower Monte Carlo in order to produce the final state observables (hadron 
level). We used PYTHIA [23] to model the parton evolution and hadronization of quarks and 
gluons emitted during the BH evaporation, as well as to treat the unstable particle decays and 
the proton remnants. This produced the final features of the events as would be seen by the LHC 
detectors. Besides the signal BH samples, also one of the most relevant SM backgrounds for such 
topologies, the process pp — ttX, has been generated using the NLO generator POWHEG [30] . 
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pp collision ^ = 7 TcV 


Process 


Mq [TeV] 


D 


a [pb] 


remnant 


P^'^P > 0.5 TeV 


P™'=" > 0.5 TeV 


Et>2 TeV 


BH 


2.0 


6 


0.017 


KINCUT=TRUE 


0.18 ±0.02 


11.8 ±0.1 


0.39 ± 0.03 


BH 


2.0 


6 




KINCUT=FALSE 


0.51 ± 0.03 


14.9 ±0.2 


7.1 ± 0.1 


BH 


2.0 


6 




no remnant 


0.94 ± 0.04 


1.74 ± 0.05 


15.8 ±0.2 


BH 


2.0 


8 


0.052 


KINCUT=TRUE 


0.49 ± 0.05 


33.9 ±0.4 


1.69 ±0.09 


BH 


2.0 


8 




KINCUT=FALSE 


1.57 ±0.09 


43.0 ±0.5 


22.9 ±0.3 


BH 


2.0 


8 




no remnant 


3.2 ±0.1 


5.4 ±0.2 


46.5 ± 0.5 


BH 


2.5 


8 


1.74 X 10-^ 


KINCUT=TRUE 


0.0023 ± 0.0002 


0.130 ±0.001 


0.0169 ± 0.0005 


BH 


2.5 


8 




KINCUT=FALSE 


0.0070 ± 0.0003 


0.151 ±0.002 


0.112 ±0.001 


BH 


2.5 


8 




no remnant 


0.0119 ±0.0005 


0.0225 ± 0.0006 


0.162 ±0.002 


tt 






84 







2.0 ± 1.2 


5.4 ±1.9 


pp collision ^ = llToV 


Process 


fie\ J 


D 


a tpbj 


reimiaiil 


7'.)^ '' > l Te\" 


I'i^."-'-' > 1 Te\' 


El > ITcV 


BH 


3.5 


6 


0.040 


KINCUT=TRUE 


0.30 ± 0.03 


24.8 ± 0.3 


0.47 ± 0.04 


BH 


3.5 


6 




KINCUT=FALSE 


1.01 ± 0.06 


32.7 ±0.4 


10.4 ±0.2 


BH 


3.5 


6 




no remnant 


1.82 ±0.08 


3.67 ± 0.12 


33.5 ±0.4 


BH 


4.0 


6 


3.6 X 10"^ 


KINCUT=TRUE 


0.037 ± 0.004 


2.45 ±0.03 


0.089 ± 0.006 


BH 


4.0 


6 




KINCUT=FALSE 


0.102 ± 0.006 


3.02 ±0.03 


1.44 ±0.02 


BH 


4.0 


6 




no remnant 


0.198 ±0.008 


0.36 ±0.01 


3.18 ±0.03 


tt 






469 










0.75 ± 0.75 



Table 2: List of generated Monte Carlo samples. Also reported are the expected numbers of events 
produced at the LHC for an integrated luminosity of lfb~^ at = 7 TeV and 14 TeV after basic 
selection cuts on the Pj- of the leading lepton (electron or muon), missing Pj- and visible transverse 
energy, for different values of Mq, and D. The expected events for one of the most relevant SM 
process, the tt production, are also reported. 



interfaced with PYTHIA for the parton shower, hadronization and particle decays. At least one of 
the two top quarks was requested to decay in either e, fj, or t. 

5 Results and outlook 

In order to point out the main characteristics of the different scenarios under study and the possi- 
bility to observe a signal over the SM background, we have applied a basic selection to the hadron 
level Monte Carlo samples using the following event variables: 

• P}^^, defined as the transverse momentum of the leading lepton (e or /x) with \r]^'^^\ < 2.5, 
where r/''^P is the lepton pseudorapidity; 

• the missing transverse momentum, 



Hp^^ +{llPy^ ' (5-1) 

where Px^ and Py^ are the cartesian x and y components of the momentum of the i^^ particle, 
and i runs over all the undetectable final state particles: neutrinos, gravitons and neutral BH 
remnant; 

• the visible transverse energy, 

= E {\fPl+Pl) ' (5-2) 
k 



jmiss 
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where k runs over all the detectable final state particles. 

The choice of the above simple variables is motivated by the fact that the observation of final state 
particles with transverse momentum of several hundreds of GeV's is the typical signal for the decay 
of a large mass state. In particular, the requirement of high-energy leptons and/or high missing 
transverse momentum allows to cope with the huge QCD background at the LHC. Since BHs decay 
democratically to all SM particles, the search for extremely energetic leptons is, in this context, 
one of the most direct way to look for deviations from the SM predictions. 

However, in case of stable neutral remnants, the lepton signal can be depressed by the fact that 
a relevant fraction of the BH mass is not available in the decay. On the other hand, the neutral 
remnant, behaving as a WIMP, will carry away a lot of energy and enhance the missing transverse 
energy signal. This is clearly visible in Figs. |3] and |4j where we show, for two representative cases at 
^/s = 7 and 14 TeV, the Pt of the leading lepton, the P^^^^ and the i?™ for the different scenarios. 
The tt background is also displayed. 

The main characteristics of the different cases are further clarified in Fig. [5j which shows, 
for a representative case at ^/s = 7 TeV, the event distribution in the E^'^-P^^^^ plane for the 
different remnant scenarios and for the tt background. If the BH decays completely, the events are 
concentrated at relatively low P™'^*^ and very high E^^^, whereas, in the case of a stable remnant, 
a shift towards much higher values of P^^^^ and consequently lower Ej^^ is clearly visible. As 
already stressed in Section [3j the control parameter KINCUT sets two extreme and probably not 
realistic descriptions of the stable remnant. Nevertheless, given the large uncertainties in the 
theoretical picture, these two borderline scenarios allow to quantify the possible spread of the 
relevant experimental observables. 

Table [2] summarizes the expected event yields, after basic selection cuts on the event variables 
described above, for the different BH scenarios and for the ti background]^ The reported values 
are normalized to an integrated luminosity of lfb~^ for y/s = 7 and 14 TeV. The results show that 
stable remnant scenarios up to Mq < 2 TeV can be already probed with the 5fb~^ of integrated 
luminosity collected by the LHC during the last year of data taking. During the year 2012, we 
expect 4-to-5 times more integrated luminosity, which will allow us to probe the fundamental mass 
in the range 2TeV< Mq < 2.5 TeV. The sensitivity to higher values could be gained only by 
increasing the LHC energy, since the production cross section for larger mass state is limited by the 
parton density functions of the proton, which steeply decrease going to larger momentum fraction 
of the partons. As shown in the second part of Table [2j Mq ~ 4 TeV could then be probed in the 
future, when the LHC runs at y/s = 14 TeV. 

Our findings should be mainly considered as an attempt to outline a reference picture, which can 
stimulate further investigations in such interesting stable BH remnant scenarios. Besides including 
specific theoretical models for BH remnants, further developments should consider other possible 
SM backgrounds (single-top, boson and di-boson production) which, even if less important than 
tt, can give non-negligible contributions. The QCD production of events with two or more jets 
should also be considered due to its huge cross section. Since such processes can contribute to the 
background only for detector effects leading to bad reconstruction of the experimental observables, 
a reliable estimation of the expected background cannot disregard a detailed detector simulation. 



''in the tt process, the high energy tails of the selection variables are not well populated due to the exponential 
decrease of the production rate with the energy. Anyway, despite possible statistical fluctuations, the contribution 
of the background in the signal region is very limited, given the very different high energy behaviour, as is clearly 
shown in Figs. [3| and [4| 
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Figure 2: Primary decay products of the BH (panels a, c, e) and multiplicity (panels b, d, f) for 
no remnant, KINCUT=FALSE and KINCUT=TRUE, respectively. The plots have been obtained 
using Monte Carlo samples generated at 7TeV, with Mq = 2TeV in D = 8. 



10 





1800„„, 2000 
P" [GeV] 




11 




10= 
1 10^ 

LU 

10 

1 

10' 
10-= 
10-^ 












[GeV] 




Figure 4: Distributions of Pr^'^, P™"^*^ and E^^^ for BH events with different remnant scenarios and 
for the ti background. The BH distributions have been obtained using the simulated samples at 
^ = 14TeV with Mq = 4TeV in L> = 6. 
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Figure 5: Distributions of signal and background events in the E^^-P^^^'^ plane for the different 
remnant scenarios. The BH distributions have been obtained using the simulated samples at 
^ = 7TeV with Mq = 2TeV in D = 8. 
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